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ABSTRACT 

We have analyzed a 90 ksec long observation of the bright nucleus of M33 
made with Chandra /ACIS. We detected low-amplitude (~ 10%) highly signif- 
icant variability on timescales of ~ 5000 sec. We also find associated spectral 
variability. The two main spectral components (a power-law with F ~ 2 and a 
multicolor disk -MCD- with kT~ 0.9 - 1.2 keV) vary in relative flux. The MCD 
temperature also increases with increasing MCD flux. The pattern of variability 
is reminiscent of (but not identical to) galactic black hole binaries. A ~ 15 Mq 
accreting black hole may explain this source. 

Subject headings: Galaxies, individuals (M33) — galaxies: nuclei — X-rays: 
individual (M33 X-8) 

February 2, 2008 

1. INTRODUCTION 

M33 is a nearby (795 kpc, van den Bergh 1991) Sc galaxy. Its nucleus is the brightest 
X-ray source in the Local Group: with a luminosity of ~ 10^^ erg/sec, it accounts for ~ 70% 
of the total luminosity of the galaxy in the 0.15-4.5 keV band (Trinchieri, Fabbiano & Peres 
1991). The nature of this source (M33-X8 in Trinchieri, Fabbiano & Peres 1991) is still 
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unclear. Because of its position, coincident with the optical center of the galaxy, it has been 
suggested that X8 could be a low-luminosity AGN (see, e.g. Markert & Rallis 1983), but 
this interpretation has two main drawbacks: no activity has been found in other energy 
bands and an upper limit of only ISOOM© has been put on the mass of the central black hole 
(Gebhardt et al. 2001). Dubus et al. (1997) have found an apparent periodicity of ~ 106 
days in a set of ROSAT data with a six years coverage, that they ascribe to a "superorbital" 
modulation in a binary system with a giant star orbiting around a black hole. However, this 
finding has not been confirmed by the subsequent observations (see e.g. Parmar et al. 2001), 
that failed to match the phase prediction. Today, the most favoured hypothesis is that M33- 
X8 is an intermediate (or even stellar) mass black hole, likely in a binary system. Recent 
studies with ASCA and BeppoSAX data (Takano et al. 1994; Parmar et al. 2001) agree in 
describing its X-ray energy spectrum with a power-law with a high energy exponential cutoff 
with ~ 2 keV e-folding energy, or with a multicolor disk black-body (MCD), both consistent 
with the spectra of Galactic black-hole candidates and of ultra luminous sources (ULXs, see. 
e.g., Makishima et al. 2000). 

In this work we examine a Chandra/ACIS-S observation of the nucleus of M33, and we 
discuss its variability and spectral characteristics. 

2. THE OBSERVATION AND DATA REDUCTION 

M33 was observed with Chandra/ACIS on July 6th 2001 (Obs Id: 2023), with an expo- 
sure time of 91.8 ksec. The observation was pointed to the star formation region NGC604, 
with the aimpoint on AClS-1, therefore the nucleus of the galaxy is ~ 12' off-axis and falls 
on AC1S-S3. Though the image of the nucleus is highly degraded, this proves to be an ad- 
vantage in terms of spectral and temporal analysis, as the photons from the nuclear source 
are not affected by pile-up effects, a major problem if the source were on-axis. 

Due to the degradation of the point spread function (PSF) with increasing off-axis 
angle, the source appears distorted and elongated with an elliptical shape. The source 
counts were thus extracted from an elliptical region centered on the position of M33-X8 as 
determined from ROSAT/HRI data by Schulman & Bregman (1995), i.e. RA=0l'*33'"50.92^ 
Dec=30°39'36.7" (J2000). The axes of the ellipse were determined examining the radial 
profiles along the two main direction of the PSF and chosing the radii that enclose 99% of 
the emission in each direction, namely 40 and 25 arcsec, with a major axis position angle 
of 53°. The background was extracted from a concentric elliptical annular region with 110 
and 70 arcsec major and minor semi-axes respectively, excluding the portions of the annulus 
falling outside the chip. 
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The data analysis was performed using different tools in the CIAO 2.1 software package, 
the imaging application DS9 v.2.1, the timing analysis software package Xronos v. 5. 18 and 
the spectral analysis software XSPEC vll.2. 

3. TEMPORAL ANALYSIS 

The observation is continuous with a time coverage of 91.8 ksec, an average count rate 
of ~ 1.9 counts/sec in the 0.2-8.0 keV band and a time resolution of 3.24 s, that allows 
for a detailed temporal analysis with moderate resolution. The light-curve of M33-X8 has 
been binned in 40 intervals, each 2250 s long. The source is evidently variable: the light 
curves show a significant variability of more than 10% from maximum to minimum (Figure 1- 
Panel A). Panels B, C, D in Figure 1 show the light-curve in three energy sub-ranges (0.2- 
1.5 keV, 1.5-4.5 keV, 4.5-8.0 keV), chosen to search for energy-dependent variability: the 
soft band samples the spectral region that hosts most of the absorption edges and spectral 
lines, the medium and hard bands cover the continuum and their ratio measures the slope 
of the continuum itself. We applied to each band (including the full band) various tests: the 
Kolmogorov-Smirnof test (see, e.g., Conover 1971), the test, the CoUura test (a phase 
averaged test plus an interpretation with a pulse model, Collura et al. 1987). These tests 
(see Table 1) confirm the presence of variability with high statistical significance in the full 
band and in the soft and medium bands. The Collura test indicates a variability timescale 
of ~ 5000 sec for the total, soft and medium bands. 

The hardness ratio plotted in Figure 1 (Panel E) is calculated as HR = where 
Hard (H) is the number of counts in the 1.5-8.0 keV band and Soft (S) is the number of 
counts in the 0.2-1.5 keV band. This curve strongly suggests that some spectral variability 
occurred during the observation. A test against the hypothesis of constant value (equal 
to the average value) gives x^/^ — 57.5/37, with a 2% probability that the hardness ratio is 
constant. As suggested by Figure 1 and, more strongly, by Figure 2, there is a correlation 
between the medium -|- hard band flux and the hardness ratio. This suggests that the 
variabihty is stronger in these bands than in the soft one. We find no evidence of correlations 
between total flux and hardness ratio (Figure 2, bottom panel), and X-ray colors with 
changing fiux (Figure 3). Colors are defined as Rl= Medium/Soft and R2=Hard/Medium, 
with Soft, Medium and Hard energy bands as in Table 1 (see also Figure 1). 

A power spectrum analysis performed on the photon counts with a 3.24 s resolution in 
the frequency range 10~^ — 0.3 Hz yields a fiat frequency distribution that can be fit with a 
single power law with a slope of ~ 4 x 10~^ with x^/^ = 64.7/75, and is consistent within 
2(7 with the Poisson noise in the frequency range sampled by the observation. 
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4. SPECTRAL ANALYSIS 

4.1. Integrated spectrum 

The source spectrum appears to be affected by a PHA to PI conversion problem^ tfiat 
causes spurious emission and absorption features between ~ 1.5 and ~ 2.5 keV. Pending 
resolution of this problem, the spectral analysis has been conducted excluding the energy 
range affected by this problem. This left us with two data ranges: 0.3-1.5 keV and 2.5-8.0 
keV, that were analyzed simultaneously. Moreover, the spectrum has been corrected for the 
quantum efficiency degradation of ACIS at low energies, using the ClAO tool CORRARF. We 
find that the total spectrum can be described with a three components model: a powerlaw 
with r ^ 2.0, a MCD with kT^^j ^ 1.0 and a Raymond-Smith model (RS) with kT-- 0.18 
and solar abundance. The residuals, however, still show an excess at ~ 1 keV, that can 
be fitted with a gaussian at ~ 0.9 keV with a ~ 0.12 keV. After this new component, the 
reduced becomes 1.27, with an improvement of 82 over 3 additional d.o.f., yielding an 
F-test probability lower than 0.1%. This bump could be the signature of an overabundance 
of Fe and/or Ne in an environment of warm matter. We thus tried to model it by allowing 
the Fe and Ne abundance to vary in the RS model, but the model was not able to reproduce 
it: the best result has been obtained having the Ne abundance four times higher than the 
solar value, but the residuals still show a narrow gaussian-like feature at 1 keV that cannot 
be accounted for by varying the abundance of other elements. Thus, the feature cannot be 
associated unambiguosly to any of the continuum components. Figure 4 shows the best fit 
model residuals, including also the range that has been excluded from the fit. It also shows 
that there is no evidence of any systematic excess in the residuals, thus the value of the 
reduced can be entirely ascribed to the statistic scattering of the data. The full set of 
best fit parameter values is reported in Table 2. The total flux is Fo.s-s.ofeey = 2.35 x 10~^^ 
erg sec~^cm~^, that, at the distance of M33, gives a luminosity of 1.5 x 10^^ erg s~^. 

The fiux of the RS component is roughly consistent with the value of the 20" extended 
emission reported by Schulman & Bregman (1995) with ROSAT/HRI. To investigate this 
possibility further, we extracted the spectrum of the circumnuclear region of M33 X-8 from 
one of the publicly archived on-axis observations of this source (Obs Id: 786, Exposure time: 
50 ksec), using the same elliptical region described in Section 2 but excluding a circle with 
radius 4.92" (lOpixels) centered on the nuclear source (heavily piled up in this observation). 
The resulting spectrum (~ 3200 counts) does not match the RS in either shape or intensity. 
It can be well described by a power-law with F = 1.4, with a fiux of 4 x 10~^^ erg sec~^cm~^in 
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the 0.3-8.0 keV band (Lo.3-8.0 = 2.5 x 10=^^ erg/s), ~l/4 of the RS flux. So, most of the RS 
component originates from a region with a radius smaUer than ~ 5", corresponding to ~ 17 
pc at the distance of M33. 

4.2. Spectral variability 

Figure 1-E suggests spectral variability. In order to investigate further this variabihty, 
we have divided the data into three phases, according to the value of the hardness ratio. We 
define the '^HR-low" state as the one where the hardness ratio (HR) is lower than -0.22, the 
HR-medium" state as the one with -0.22<HR<-0.18 and the ^^HR-high" state as the one 
where HR>-0.18. We have selected three spectra according to this criterion and we have 
analyzed them separately, using the same spectral component as in Section 4.1. The results 
are in Table 2. We find that the Gaussian component is not needed for the HR-low and HR- 
high states. This is not surprising, as each of these spectra includes only ~ 10% of the total 
counts of the source, and the gaussian component is relatively faint. Examining the other 
two spectral components we find that: a) the power-law photon index F does not change 
significantly in the three spectra, while the fiux decreases from the HR-low to the HR-medium 
and HR-high states; b) the errors of the MOD components in the HR-low state are rather 
large, however there is evidence of the disk kT increasing from the medium to the hard state; 
c) we do not detect any variability in the RS, consistent with possibly extended emission 
from hot ISM. Figure 5 shows the three spectra, deconvolved from the instrumental response, 
with the best fit model components. From this figure we see a progression of larger MCD 
contribution (relatively to the power-law) going from the HR-low to the HR-high spectral 
states. 



5. DISCUSSION 

The light curve of M33-X8 shows a significant low-amplitude variability, with a ~ 10% 
variation on a timescale of a few thousands of seconds. This timescale is particularly evident 
both in the full band light curve, where the source count rate drops and suddenly goes back 
to the original value in ~ 3 x lO^s, and in the medium band light curve, where the count rate 
oscillates a few times around its average value of ~ 0.71 c/s. This result strongly confirms 
the variability on a 3000 s timescale found by Peres ct al. (1989) in the Einstein data of the 
nucleus of M33. However, the power spectrum is flat, and does not evidence any periodical or 
quasi-periodical oscillation. The hardness ratio is also variable, showing evidence of spectral 
variability, that appears to be correlated mostly with the emission at energy > 1 keV. 
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The spectral analysis reveals the presence of four components. Two of these, a Raymond- 
Smith component and a broad "Gaussian" line feature contribute relatively little to the total 
spectrum, and do not present any compelling evidence of variability. We will not discuss them 
any further. The two main components are a strong F ~ 2 power-law and a kT= 0.9 — 1.3 
keV multicolor disk component (MOD). The parameters of the power-law and of the MCD 
are in good agreement with the results found by Parmar et al. (2001), Takano et al. (1994), 
Dubus & Rutledge (2002). However, their contribution to the total spectrum is variable, 
as showed by the HR-selected spectra (Table 2): the HR-low phase is characterized by a 
brighter power-low and a fainter disk. While we do not find any further variation of the 
power-law flux from the HR-medium to the HR-high state, the MCD flux seems to increase 
with increasing kT. This result may suggest that the variability is linked with fluctuations 
in the accretion rate. 

For the MCD component, the lower limit of 37 km to the inner radius of the disk 
(Table 2; the limit is due to the presence of the V cos 6 factor accounting for the inclination 
angle of the disk) implies a lower limit to the mass of the compact object of ~ 4.2Mq. The 
Eddington limit for this mass is LEdd = 6.3 x 10^^ erg/sec, that is consistent with the power- 
law + disk emission of the source. This is in agreement with the conclusion of Parmar et 
al. (2001), Takano et al. (1994) and Dubus & Rutledge (2002) that M33 X-8 is most hkely a 
normal binary system containing a stellar mass black hole that happens to be at the center 
of the galaxy, rather than an AGN. Although not extreme, the spectral variabihty of M33 
X-8 is reminiscent of what has been observed in black hole binaries. In the Galactic black 
hole binaries we usually observe a "high/soft state", dominated by a hot disk, sometimes 
associated with a faint power-law tail, and a "low/hard state" dominated by a power-law 
spectrum, with a colder disk component (see, e.g., Esin, Mc Clintock and Narayan 1997). 
The soft state is explained with high accretion rate, with the accreting disk going down 
toward the last stable orbit, while the hard state is associated to a lower accretion rate that 
causes the disruption of the inner part of the disk and the presence of a comptonized corona 
wrapping the central source. In the case of M33 X-8 both disk and power-law component 
are always present, but the power-law dominated HR-low state corresponds with a lower flux 
in the medium energy band (1.5-4.5 keV). With increasing medium band flux the spectrum 
evolves to be dominated by the disk component. 

The source that resembles more M33 X-8 is probably LMC X-3 (Nowak et al. 2001), 
a black hole candidate with a probable mass of 9Mq and a luminosity of up to 4 x 10^^ 
erg/s. M33 X-8 presents many similarities with LMC X-3: the thermal spectrum, that can 
be modelled with a MCD (with kT~ 0.8 in LMC X-3 and kT~ 1.2 in M33 X-8) and a steep 
power-law; the presence of variability on timescales of some ten thousands of seconds; a long 
timescale periodicity (116 d for M33 X-8 and 99 or 198 d for LMC X-3, see Dubus et al. 
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1997 and Cowley et al. 1991) that appears to vary or even disappear in a few years (Nowak 
et al. 2001; Parmar et al. 2001). Also, LMC X-3 is characterized by a lack of fast variability 
(on timescales < 1 s), but this feature cannot be investigated in M33 X-8, due to the lack of 
data with high enough temporal resolution. However the power-law component is relatively 
more prominent in M33 X-8. 

Concluding, the Chandra spectral variability observations of the nucleus of M33 (X-8) 
strongly agree with the suggestion that this source is an accreting > 5 Mq black hole with 
a variable accretion disk, comparable to Galactic black hole binaries. 

We thank M. Krauss and F. Nicastro for their help in solving the calibration problems, 
J. MacDowell for interesting discussions and A. Kong, D. Schwartz and J. Drake for useful 
comments on the manuscript. This research has been supported by CNAA and by NASA 
contract NAS 8-39073 (CXC). 
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Fig. 1. — Background subtracted light curves of M33-X8. Time is from the beginning of the 
observation. Panel A: total Chandra hand (0.2-8.0 kcV). Panel B: light curve in the 0.2-1.5 
keV energy band. Panel C: light curve in the 1.5-4.5 kcV energy band ; Panel D: light curve 
in the 4.5-8.0 keV energy band. In all panels the average value is marked as a solid line. 
Panel E: hardness ratio calculated using the 1.5-8.0 keV and 0.2-1.5 keV energy bands. Bins 
are 2250 s long. 



-10- 




Fig. 2. — Count rate in the two bands (S and H, used for the hardness ratio) vs. hardness 
ratio and total count rate vs. hardness ratio 
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Fig. 3. — Color-color diagram (the three energy bands are defined in Table 1. The filled 
points are those occurring during the low fiux regime (i.e. the dip between 5.0 x 10^ s and 
7.5 X 10^ s; see Figure 1). 



Table 1 . Variability of the nucleus of M33 





Full band 


Soft band 


Medium band 


Hard band 




0.2-8.0 keV 


0.2-1.5 keV 


1.5-4.5 keV 


4.5-8.0 keV 


Avg count rate 


1.96 ±0.05 c/s 


1.15 ±0.04 c/s 


0.72 ±0.03 c/s 


0.052 ±0.008 c/s 


RMS frac.var.(«) 


2.5 % 


2.2% 


3.3% 






148.8/37 


98.8/37 


107.8/37 


39.1/37 


KS probability^'^) 


1.1 X 10-5 


5.9 X 10-6 


1.3 X 10-3 


0.06 


Timescale^*) 


- 5000 s 


~ 5000 s 


~ 5000 s 





'^Root Mean square fractional variation (only where higher than the 3a level) 

''calculated against the the hypothesis of having a constant count rate equal to the 
average value. 

'^Kolmogorov-Smirnov test probability against the hypothesis of constant rate. 
^Variability timescale from the Collura test (CoUura et al. 1987) 
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Fig. 4. — Spectrum of M33-X8 and residuals over the best fit model (Table 2). The energy 
channels between the two vertical lines (1.5-2.5 keV) have been excluded from the fitting 
because of an abound conversion problem (see Section 2) 
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KeV 

Fig. 5. — Data and best fit models (unfolded for the instrumental response) for the 
three phase-selected spectra, deconvolved from the instrumental response. PL=Power-law; 
MCD=Multicolor black body disk; RS=Raymond-Smith; G=Gaussian 
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Table 2. Spectral Analysis - Best fit parameters 



Component 


Parameter 


Integrated emission^^^ 






Total spectrum HR-Low phase Medium phase HR-High phase 





Nh (xl02icm-2) 








9 1+0.4 


2 0+0-2 


2.1 ±0.6 


Power law 


r 






2.18 ±+0.07 


2 14+011 


9 1 9+0.11 
■^•1-^-0.07 


9 c+0.2 

^•"3-0.3 




Flux (10~^^ erg sec" 


^cm 




1.45 ± 0.09 


1.9 ±0.2 


1.4 ±0.1 


1.3 ±0.2 


MCD 


kTi„t (keV) 






1.14 ±0.05 


1 o+°-3 

■'■•^-0.2 


1 1 1 +0.06 
-•-■11-0.07 


1.29 ±0.08 




R\/cos6' (km) 






391^ 




42 ±4 






Flux (10~^^ erg sec" 


^cm" 




8.5 ± 1.8 




9±2 




Raymond-Smith 


kT (kcV) 






n 18+0-01 


0.19tn^ 


0.17 ±0.02 


0.20 ±0.05 




Flux (10^^^ erg sec" 


■^cm~ 




1.7 ±0.4 


1 7+1-1 


1.8 ±0.9 


1 l+'i-^ 
-'-•I--0.7 


Gaussian 


E (keV) 
cr (keV) 






n in+0 05 




04+0.04 

^•^+-0.10 

n 11+0-09 

^•11-0.04 






Flux (10"^^ erg sec" 


^cm~ 




1 7+4.3 
^- ' -0.6 




9 9+3.3 
^•^-0.6 












434/333 


210/203 


363/314 


182/198 



Circumnuclcar emission^^^ 
(xlO^icm-^) 0.6 ±0.3 

Power law F 1.40 ±0.10 

Flux (10-13 erg sec-icm"^) 4.7 ±0.5 
X^/i^ 54/66 



(1) Prom Obs Id 2023 

(2) Prom Obs Id 786; see Section 4.1 



Note. — Unabsorbcd fluxes arc calculated in the 0.3-8.0 energy band. 9 is the inclination angle of the disk in the MCD 
model. The Gaussian component was not needed to fit the low and high phases. Quoted errors are at 90% confidence level 
for one interesting parameter. 



